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1 Problem 

This note considers variations on the theme of a solenoid magnet {i.e., a magnet whose field 
has axial symmetry) as a lens for charged particles. A related problem has been posed in 

III 

Recall that if a device is to be a lens with optic axis along the z axis in a cylindrical 
coordinate system (r, 0, z) , then as particles leave the device they must have no azimuthal 
momentum, P^ = 0, and their radial momentum must be proportional to their radial coor- 
dinate, Pr oc r. Special cases are (1) that all particles have r = at the exit of the device, 
which is a focal point; and (2) that all particles have zero radial momentum. 

1.1 Particle Source Inside the Solenoid: A Neutrino "Horn" 

A neutrino "horn" is a magnetic device whose goal is to focus charged vr mesons that emerge 
from a target into a parallel beam, so that when the pions decay, vr^ — >■ /i^z/, the resulting 
neutrinos form a beam that has minimal angular divergence.^ Suppose the pions are pro- 
duced at the origin, inside a solenoid magnet of uniform field B = Bz whose axis is the z 
axis and whose downstream face is at 2; = L. Show that pions of momenta 

eBL , , 

fn = 0,l,2 ) (1) 



{2n + l)7rc 



emerge from the magnet with their momenta parallel to the z axis, independent of the 
production angle 9 (for 9 <ti 1). In this case, the solenoid acts like an ideal thin lens of focal 
length L, located at z = L. 

Neutrinos from the forward decay of the resulting parallel beam of pions will have a quasi 
line spectrum with momenta proportional to those of eq. (^. If the neutrinos are detected 
at a distance / from the source, that distance can be chosen so that the various peaks in 
the neutrino spectrum all satisfy the condition for maximal probability of oscillation into 
another neutrino species prior to their detection. 

1.2 Particle Source Outside the Magnet 

Consider a point source of charged particles located at a distance D from the entrance to 
solenoid magnet of length L and field strength B, the source being on the magnetic axis. For 
what momenta P are particles with angle ^ <^ 1 with respect to the magnetic axis focused 
to a point on axis beyond the exit of the magnet? 



^Because of the Jacobean peak in the two-body decay kinematics of the pion, for some purposes it is 
favorable to use neutrinos produced at a nonzero decay angle. See, for example, 
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In both cases, the focusing effect is due to the fringe field of the magnet, and not due to 
the uniform central field. A simple model of this effect (impulse approximation) supposes 
the magnetic "kicks" of the fringe field occur entirely in the entrance and exit planes of the 
magnet. Although this effect can be analyzed by direct use oi F = ma, it is helpful to 
consider the canonical (angular) momentum of the particle in the magnetic field. For this, 
you can use either a Lagrangian formulation, or direction calculation via the Lorentz force 
law, in which latter case first consider dL^/dt = d{r x P)z/dt. 



2 Solution 

Although this problem can be solved without explicit use of the canonical angular momen- 
tum of a charged particle in a magnetic field, that concept offers an elegant perspective. 
Therefore, we first discuss canonical momenta in sec. 2.1, and then comment on the paraxial 
approximation in sec. 2.2, and the impulse approximation in sec. 2.3, before turning to the 
solutions for solenoid focusing of particles produced outside, and inside, of the magnet in 
sees. 2.4 and 2.5. The possibly novel aspect of this note is the discussion in sec. 2.5.1 of a 
neutrino horn based on solenoid focusing. 



2.1 Conservation of Canonical Angular Momentum 

The canonical momentum of a particle of charge e and rest mass m is (in rectangular coor- 
dinates and in Gaussian units) 

P = P + — , (2) 



where P = 7mv = mv/y 1 — v'^/c^ is the mechanical momentum of the particle, A is the 
vector potential of the magnetic field, and c is the speed of light. The canonical angular 
momentum is 

1 = r X p, (3) 

where r is the position vector of the particle. 

One way to deduce the conserved quantities for the particle's motion is to consider its 
Lagrangian or Hamiltonian. If an electric field is present as well, with electric potential V, 
the Lagrangian C of the particle can be written [3] 

^ mc^ e A ■ V ^ ^ , ^ , 

£ = + eV, (4) 

7 c 

where v = dr/dt is the particle's velocity. The canonical momentum associated with a 
rectangular coordinate Xi is therefore pi = dC/dxi, leading to eq. (j2I). Then, the Hamiltonian 
Ti of the system is 



n 



\ 



A \ 2 

eA 



m^c'' + p ] +eV. (5) 



If the external electromagnetic fields have azimuthal symmetry, then the potentials V 
and A do also. We consider a cylindrical coordinate system (r, 0, z) with the z axis being 
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the axis of symmetry of the fields. Then both the Lagrangian and the Hamiltonian have no 
azimuthal dependence, 

'4-^ = 0, (6) 



SO the equations of motion (and the identities r = rf + zz, r = v = rf + r(f)(f) + iz) tell us 
that the canonical momentum is a constant of the motion (even for time-dependent fields, 
so long as they are azimuthally symmetric),^ 



dC 



r ( ^mrcp + 



L. 



(7) 



We also see that the canonical momentum can be interpreted as the z component of the 
canonical angular momentum (jS)), so Iz is also a constant of the motion. 

For completeness, we verify that dlz/dt = using the Lorentz force law. 



— = e(E+-xB 

at 
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c ot c 



We begin with the ordinary angular momentum L = r x P, and consider the z component 
of its time derivative: 



dLz _ d{r X P) 

~dr ~ 



dt 



r X 



dP\ 



'rfP\ 



(9) 



From eq. 

Vdl), 



we have, since dV/dcj) = dAr/dcj) = dAz/dcj) = 0, 
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(10) 



where ^ when applied to a field such as the vector potential A is the convective derivative 
associated with the moving particle. 

Noting that P = 7m(rr + r(j)(f) + zz) and r = 00, we also find 



'dP\ 
It) 



dP^ 
dt 



+ 



Combining eqs. ©-(HD), we have 

dL, 



Hence, 
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dt 
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dlz 
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dt 



(12) 
(13) 



as found by the Lagrangian method as well. 

^Note that the definition Q of the canonical momentum leads to the awkward result that — r{p)^, 
where (p)^ is the (j> component of the canonical momentum vector p of eq. 1^)1. 
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2.2 The Paraxial Approximation 

We now turn our attention to the question of lenslike character of a solenoid magnet as a 
charged particle moves from a region of uniform field to zero field, or vice versa. 

Inside a uniform solenoidal magnetic field B = Bz, the trajectory of the particle is a helix 
(whose axis is in general at some nonzero radius tq from the magnetic axis). The radius R 
of the helix can be obtained from F = Ma = ev/c x B using the relativistic mass M = 7m. 
The projection of the motion onto a plane perpendicular to the magnetic axis is a circle of 
radius R and the projected velocity is v±. Hence, 

^ = e^5, (14) 

so that 

(15) 

where Pi = 'ymv± is the transverse momentum of the particle. For a particle whose average 
velocity is in the +z direction, the sense of rotation around the helix is in the —(f) direc- 
tion (Lenz' law). The angular frequency u of the rotation (called the Larmor or cyclotron 
frequency) also follows from eq. ()14p: 

u = — = . (16) 

R 'ymc 

If the solenoid magnet has length L, then the time t required for the particle to traverse 
the magnet is given by 

* = - = ^ = T^. (17) 

where 6 is the production angle of the particle with respect to the z axis. Hence, the 
trajectory of the particle rotates about the axis of the helix by azimuthal angle (ph as the 
particle traverses the magnet, where 

eB jmL eBL 

^h = uJt= ^ - = — -. 18 

7mc R cos o cR cos o 

There is a unique value for (ph only for small production angles {6 <^ 1), which is called the 
paraxial regime: 

eB L 

(ph ~ — L = — , (paraxial approximation, 9 <ti 1), (19) 
cP X 

where we define the (reduced) Larmor wavelength of the particle's motion to be 

cP 

A.-. (20) 

In the paraxial approximation the magnetic force that bends the particle's trajectory into 
a helix is a weak effect, in that it depends on the product of the small transverse velocity 
v± = V sin 6 <^v and the axial field B. 
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2.3 The Impulse Approximation 

As the trajectory crosses the fringe field of the solenoid, the axial field drops rapidly from B 
to zero (or rises rapidly from zero to B). In this region there must be a radial component 
to the magnetic field, according to the Maxwell equation 



ld{rBr) , dB, 
r r dz ^ 



= V-B = -^^^l^ + ^, (21) 



so that 

-—it <-) 

(as also readily deduced by applying Gauss' law to a "pillbox" of radius r and thickness dz). 
Although the radial component B^. of the magnetic field is small, it couples to the large axial 
velocity to give a force = dP^/dt in the azimuthal direction that is not negligible. We 
can write 

dP^ 1 dP^ 1 eVzBr erdB^ 
dz Vz dt Vz c 2c oz 

Hence, the change AP<^ in the azimuthal momentum of the particle as it crosses the fringe 
field is 

erABz erAB 

APs ^ = , 24 

2c 2c ' ^ ^ 

since ABz = —B at the axial field falls from B to zero. 

The impulse approximation is that during the particle's passage through the fringe 
field we can neglect the change in its momentum due to coupling with the axial magnetic 
field. We only consider the azimuthal kick (j24j) . Thus 

er B 

Pr,out = -Pr.in, -P</.,out = P<f,;m + , -P^.out = -Pj.in (impulse approximation). 

(25) 

Furthermore, we neglect the change in the transverse coordinates of the particle as it passes 
through the fringe field. 

?"out = '"in, 0out = 0in (impulse approximation). (26) 

We can connect the impulse approximation with conservation of canonical angular mo- 
mentum by noting that a solenoid magnet with (uniform) field B = i?z has vector potential 

A = 71^0 = ^0. (27) 

To see this, recall that B = V x A implies that the integral of the vector potential around 
a loop is equal to the magnetic fiux through the loop; hence, 27rrA^ = nr^B. 

The z component of the canonical angular momentum (which is equal to the azimuthal 
component of the canonical momentum 



lz=P4> = r{P^ + eA^/c) = r{P^ + erB/2c), (28) 
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is a constant of the motion for a particle in a solenoid magnet. Hence, we see that the 
simplified impulse approximation that rout = '^in plus conservation of canonical angular 
momentum implies the form (|25p. 

Additionally, we note that particles which are created on the magnetic axis have Iz = 0, 
whether they are created inside or outside the magnetic field. As a consequence, whenever 
such a particle is outside the magnetic field region it has = 0. If it has passed through 
a region of solenoidal magnetic field, the azimuthal kicks at the entrance and exit cancel 
exactly. This results does not depend on the impulse approximation, as it is deduced directly 
from conservation of canonical angular momentum. 



2.4 Particle Source Outside the Magnet 

We consider a solenoid magnet whose axis is the z axis with field B = Bz for < 2 < L. 
A particle of momentum P and charge e is emitted at polar angle ^1 <C 1 from a (point) 
source at {x,y,z) = {0,0, —di), and so arrives at the entrance of the magnet with spatial 
coordinates (r, (p, z) ~ (ri = diOi, 0, 0) in the small angle (paraxial) approximation, and with 
momentum [Pr, P^, Pz) ~ {Pr^,0,P), where 

Pr, = P^l. (29) 

The projection of the particle's trajectory onto the x-y plane is shown in Fig. ^ 

The fringe field at the entrance of the solenoid gives the particle an azimuthal kick 
resulting in momentum 

P,. = -^ = -^^, (30) 
2c 2c ' ^ ' 

according to eq. where the magnetic field is B = Bz inside the solenoid. The transverse 
momentum Pi of the particle inside the magnet is therefore 



eBri 



f 2cP \ ^ eBri 



2^^' = ^. (31) 



where R is the radius of the helical trajectory of the particle inside the solenoid, recalling 
eq. psp . We also can write 

ri = 2Pcosa, (32) 
where the angle a, shown in Fig. ^ is related by 

-Pri 2cP 2A 
tana = ^— = — — = — , (33) 
IP^J eBdi di 

which is independent of the production angle 61 in the paraxial approximation. 

As the particle traverses length L of the solenoid, its trajectory rotates by azimuthal 
angle 

eBL L , , 

about the axis of the helix. At the exit of the solenoid the particle is at {ji, (f), L) in cylindrical 
coordinates centered on the axis of the magnet (rather than on the axis of the helix), as shown 
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Magnetic Axis 




Figure 1: Geometry of the helical trajectory of a particle of total momentum 
P that enters a solenoid magnet at {r,(j),z) = (ri = di6i,0,0) with radial 
momentum = P6i. The fringe field at the entrance of the solenoid gives 
the particle an azimuthal kick resulting in momentum P^^ = —eBri/2c, where 
the magnetic field is B = Bz inside the solenoid. The helix has radius R = 
cP±/eB. At the exit of the solenoid the particle is at (r2,0, L) where = 
—eBL/2cP = 0/i/2; the azimuthal rotation of the particle's trajectory about 
the magnetic axis is one half that about the axis of the helix. 



in Fig. nj By the well-known geometrical relation that the angle subtended by an arc on a 
circle as viewed from another point on that circle is one half the angle subtended by that 
arc from the center of the circle, we have that^ 

(ph eBL L 

The radial coordinate of the particle at the exit of the solenoid is 

r2 = 2Rcosp, (36) 



•^The geometrical relation ()35|1 has the consequence that in a frame that rotates about the magnetic 
axis at half the Larmor frequency Hlt)|) . the particle's trajectory is simple harmonic motion in a plane that 
contains the magnetic axis However, we do not pursue this insight here. 
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where angle (3 is given by 



[5=\<P\-a = ^-i^n (37) 




When the particle is at the exit of the solenoid, but still inside it, the transverse mo- 
mentum vector P_|_2 makes angle (5 to the unit vector 0, as shown in Fig. ^ The radial 
momentum of the particle P^^ at the exit of the magnet is therefore 

using eqs. and (jHSJ, while the azimuthal component obeys 

= -P± cos/? = . (39) 

As the particle exits the magnet, the radial component of its transverse momentum 
remains at the value of eq. (jHH|l in the impulse approximation, while the azimuthal component 
increases by ePr2/2c over the value of eq. and hence vanishes, as expected since the 
canonical angular momentum is zero. 

Once the particle has exited the magnet its transverse momentum is purely radial, with 
a value proportional to the radial coordinate r2 at the exit of the magnet. This is lens-like 
behavior, in that the particle will then cross the magnetic axis at distance ^2 from the exit 
of the magnet, where 



and so 



where 



rf. = ^H^ = 2Ai±l^ = ^fl + H^ta„A|, (41) 

eP tan/3 tan ^ - f di - f \ di 2X' ' ^ ' 



tan ^ 

When distance d2 is positive the solenoid acts as a (thick) focusing lens. 

For the special cases of point-to-parallel focusing (^2 oo) and parallel-to-point focusing 
{di oo), the solenoid magnet has focal length / given by eq. 

If {2X/di) tan(L/2A) ^ 1 then the object distance di and the image distance d2 obey the 
lens formula 

1 1 1 / L di\ 

* + 4 = 7 r"5A«23;j^ 

If in addition the length L of the solenoid is small compared to the Larmor wavelength 
A the solenoid can be called a thin lens, for which 

f = — (thin lens : L < A, L < di). (44) 
ij 

This weakly focusing limit is, however, seldom achieved in practical applications of solenoid 
magnets as focusing elements. 

The results (jlH)-(03) thick- and thin- lens focusing can be utilized in a transfer-matrix 
description of particle transport through magnetic systems [H]. 
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2.5 Particle Source Inside the Magnet 



The case of a source of particles inside the solenoid magnet, say at 2; = 0, can be treated as 
a special case of the analysis in sec. 2.4 in which di = ri = 0. The angle a shown in Fig. ^ 
is 7r/2 in this case, so that angle P is 

The radial coordinate of the particle at the exit of the magnet is 

r2 = 2i?cos/3 = 2i?sin;^ , (46) 

2A 

and the image distance d2 follows from eq. ()4H) as 

2cP ^ L 

d2 = — ^ = -2Atan — . (47) 

eBtan(3 2A ^ ' 

The radial momentum at the exit of the magnet is 

^ „ ^ eBR To ^ eBro L , , , 

p = -p^ sin P = ^ tan /3 = cot — . 48 

' ^ c 2P 2c 2A ^ ' 

according to eqs. and 

This is lens-like behavior (Pr2 oc r2) for any length L of the solenoid, with L = nnX being 
the boundary between focusing and defocusing. 

For the special case that L = 2miX we have ^2 = ^"2 = 0, corresponding to an image of 
the source occuring at the exit of the magnet. 

Of particular interest here is the special case that L = {2n + l)vrA, for which d^ = 00, 
Pr^ = 0, and we have point-to-parallel focusing. From Fig. [T] and eq. (jl^ we see that the 
condition for point-to-parallel focusing of a source inside the solenoid is that the particle has 
completed an odd number of half turns on its helical trajectory when it reaches the end of 
the solenoid. In this case we can say that the focal length of the solenoid lens is just the 
length L, 

cP 

f = L = {2n + 1)ti—- (point-to-parallel focus, source inside solenoid). (49) 
eB 



2.5.1 Neutrino Horn: Point-to-Parallel Focus, L = (2n + l)ncP/eB 

A solenoid magnet provides point-to-parallel focusing for particles produced inside the mag- 
net, on its axis, with a discrete set of momenta P„ given by 

^" = Tr^' (n = 0,l,2,....) where Pq = ^ . (50) 
2n + 1 ire 

Particles with other momenta are not brought into parallelism, so that a "beam" formed 
by drifting particles that emerge from the solenoid will be quasimonochromatic with the 
sequence of momenta given in eq. (j^OI)- Figure |21 illustrates trajectories for particles of 
momenta Pq and 3Po in a solenoid magnet. 



9 




Figure 2: Concept of a neutrino horn based on solenoid focusing. The pion 
production target is inside the uniform field region of the solenoid. The fo- 
cusing effects of the fringe field at the exit of the magnet (at distance L from 
the target) act as ideal thin lens of focal length L for a discrete set of particle 
momenta, given in eq. ()50|) . 



Such a sequence of momenta occurs in the phenomenon of neutrino oscillations over a 
flight path /. As is well, known, in the approximation of pure two-neutrino mixing, the 
probability that neutrino type (mass eigenstate) i of energy E = P appears are neutrino 
type j after traversing distance / is given by 

Prob(i j) oc sin^ , (51) 

2E 

where AMij = Mi — Mj is the difference in the masses of the two neutrino types. Hence, for 
a fixed drift distance /, the probability of neutrino type i appearing as type j is maximal for 
the sequence of neutrino momenta (energy) 

P AM ^ / 

^" = Tr^' (^ = 0,1,2,....) where Pq = ^- (52) 

2n + 1 TT 

Thus a solenoid magnet could be very useful in preparing a neutrino beam with a sequence 
of momenta such that all oscillation effects are maximal. The potential advantage of such a 
beam for the study of CP violation in neutrino oscillations has been pointed out by Marciano 
[H], and elaborated upon in [7j. 

Of course, neutrinos are neutral, so that a solenoid magnet cannot directly affect their 
trajectories. Rather, the solenoid magnet would be used to focus vr^ particles that are 
produced in the interaction of a proton beam with a nuclear target that is placed on the axis 
inside the magnet. The length / of the magnet should be short enough that most pions of 
interest exit the magnet before decaying into neutrinos, according to 

7r+ li^i^^i, tt" /i"z/^. (53) 

Because of the low "Q" value of this decay, the direction of the neutrinos is very close to 
that of the pions, provided that latter have energies greater than a few hundred MeV. The 
forward-going neutrinos carry about 4/9 of their parent pion momentum, so the solenoid 
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system should be chosen with a momentum Pq.tt equal to 9/4 of the highest desired neutrino 
momentum at which the oscillation probability is maximal, i.e., 

Po,. ~ Ipo,.. (54) 

As implied by eq. the solenoid-focused beam would contain both muon neutrinos 
and muon antineutrinos, in roughly equal numbers. This has the advantage to studies could 
be made simultaneously with both neutrino and antineutrino beams. However, for the study 
of CP violation it would be necessary to identify whether each interactions was due to a 
neutrino or an antineutrino. This identification must be provided by the detector in which 
the neutrino interacts. If the neutrinos oscillate into electron neutrinos or antineutrinos 
before they interact in a the detector, the latter must distinguish showers of electrons from 
positrons. This difficult experimental challenge can likely only be met by a magnetized liquid 
argon detector [HI El CHI • 

When studying the oscillation of muon neutrinos into electron neutrinos, the presence of 
electron neutrinos in the beam constitutes the limiting background. Electron neutrinos are 
present in the beam due to the 3-body decay of the muons from pion decay: 

7r+ fi^i^t,, e^UeUf,, IT' f^'i^ti, /i~ e'Uei^f,. (55) 

The background of electron neutrinos, compared to the flux of muon neutrinos at a particular 
energy, is suppressed when the beam contains only a narrow range of momenta of the parent 
pions. This occurs because the muon neutrinos from the pion decay then have typically 
higher momentum that the electron neutrinos from the related muon decay. Hence, the 
solenoid- focused neutrino beam, with its quasi line spectrum of energies will have lower 
electron neutrino content, at least for highest-energy neutrino "lines", compared to a wide- 
band neutrino beam. 

A flnal advantage of the solenoid-focused beam is that the magnetic elements are farther 
removed transversely from the pion production target, and so can be made more radiation 
resistant to intense proton fluxes than is the case for more conventional toroid-focused neu- 
trino "horns". Further, the relatively open geometry of the solenoid lens will permit use of 
liquid metal target, as needed if the proton beam has several megawatts of power 

The author thanks Ron Davidson for the demonstration that conservation of the canonical 
momentum follows from the Lorentz force law. 
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